Realizing the strongly correlated rf-Mott state in a fermionic cold atom optical lattice 
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We show that a new state of matter, the d-Mott state (introduced recently by H. Yao, W. F. 
Tsai, and S. A. Kivelson, Phys. Rev. B 76, 161104 (2007)), which is characterized by a non-zero 
expectation value of a local plaquette operator embedded in an insulating state, can be engineered 
using ultra-cold atomic fermions in two-dimensional double- well optical lattices. We characterize and 
analyze the parameter regime where the d-Mott state is stable. We predict the testable signatures 
of the state in the time-of-flight measurements. 
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Experimental realizations of degenerate Bose P, Q and 
Fermi gases ^ ^, ^, JJ in cold atom optical lattices has 
led to increasing interaction between the atomic and con- 
densed matter physics communities. In fact, these atomic 
systems have become a proverbial "playground" where, 
in principle, custom Hamiltonians [J S Q can be made 
to order, enabling optical lattice emulations of strongly 
correlated condensed matter phenomena. Particularly 
exciting, in this context, is the creation of a fermionic 
Hubbard model in ultra-cold atomic systems, since the 
Hubbard model is a paradigm for theoretical studies of 
strong correlations. The two-dimensional fermionic Hub- 
bard model, a model woefully difficult to "solve", is of 
great interest since it is thought to hold the key to un- 
derstanding high-temperature (high-Tc) superconductiv- 
ity In fact, the cold atom systems allow for versatile 
tuning of the Hamiltonian parameters where the hopping 
(or tunneling) matrix element between nearest neighbor 
sites (through laser intensity tuning), the on-site inter- 
action between particles (from attractive to repulsive by 
tuning the Feshbach resonance), and the dimensionality 
(from one to three) of the model can all be varied and 
controlled Q. Such versatility in tuning the Hamilto- 
nian makes cold atom systems particularly attractive in 
studying novel quantum phases and transitions between 
them. 

The possible d-wave character of the ground state in 
the Hubbard model is an important concept in relation 
to the high-Tc cuprate superconductors |10|. Studying 
the Hubbard model experimentally in the cold atomic 
gases could, in principle, be an effective and accurate 
analog simulation of high-Tc superconductivity. Eventu- 
ally it would be very important to test whether or not 
the fermionic Hubbard modeLin its full generality, al- 
lows for d-wave superfluidity [HI, [13, [l^ , a controversial 
conjecture not yet settled theoretically. In this work, we 
study the model in a controlled limit in which there is 
(i) an exactly solvable point where a true new state of 
matter, the so-called d-Mott state recently discovered by 
Yao, et al. [l3| exists as the ground state of the checker- 
board Hubbard model (see below), and (ii) a parameter 
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FIG. 1: (Color online) (a) Two-dimensional checkerboard 
Hubbard model. Solid (dashed) bonds represent hopping 
amplitudes t {t'). (b) "Ladder" cluster we exactly diagonal- 
ize with periodic boundary conditions imposed in the "y "- 
direction and open boundary conditions in the "a;" -direction. 
The numbers label the sites of a single plaquette. 



which can be varied to destroy the d-Mott state. The d- 
Mott state |l4| is a special Mott insulating state with a 
non-zero expectation value (—1) of a local plaquette op- 
erator. Dp (see below), which cyclically rotates the sites 
by an angle 7r/2. In other words, the state is an insula- 
tor with a local plaquette d-wave symmetry-in analogy 
with the d-wave superconductor-characterized by a non- 
zero expectation value (Dp). Realizing this state in an 
optical lattice has its own intrinsic appeal because it is 
a true new state of matter not adiabatically connected 
to any other known insulating states [l^. Furthermore, 
since this is the ground state arising out of an exact solu- 
tion of the Hubbard model in a controlled limit, it would 
serve as an important benchmark for the accuracy of ex- 
periments aimed at studying the general Hubbard model. 
For this purpose, we propose and analyze how the dis- 
tinctive signatures of the d-Mott state could be observed 
in the usual time-of-ffight measurements. An alternative 
proposal was provided in Ref . [3 • 

We consider a system of fermionic atoms {'^^K) 
where the two atomic internal hyperfine states, 
\F = 9/2, mp ^ -9/2) and \F = 9/2, mp = -7/2), are 
taken as the effective spin states a ^t^i- The atoms 
are loaded into a two dimensional superlattice, produced 
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by superimposing a long and a short period lattice in 
the X- and y-direction such that an array of square pla- 
quettes is created [l^ The dynamics of the atoms 
can be described by the so-called checkerboard Hubbard 
model [17[ defined as 



H 



(1) 



Here, c^a (cj^r) is the particle destruction (creation) op- 
erator at lattice site r with spin (or magnetic sublevel) 
cr, tr{5) is the hopping amplitude for a particle at site r 
hopping to site r + 5, where site r -I- 5 is a nearest neigh- 
bor site, U is the usual on-site interaction energy, and 
Tira = cJg.CrCT is the number operator. The checkerboard 
Hubbard model is defined by choosing tr{5) = t when 5 
connects sites within a square plaquette, and tr{5) = t' 
when 5 connects sites between two neighboring square 
plaquettes, see Fig. [TJ 

When t' = the ground state of this model can be 
obtained exactly for all U /t and all band-fillings or den- 
sities (dopings) by solving the single plaquette problem 
exactly (either by brute force [3], Bethe ansatz [l9| - 
four site plaquette is a one-dimensional ring with four 
sites - or through numerical exact diagonalization) and 
constructing the full many-plaquette state as a product 
state of the plaquette basis states. At exactly half fill- 
ing (or zero doping, the case considered throughout this 
work), the ground state for U/t > is an insulator with 
(i-wave symmetry on the plaquette, the so-called d-Mott 
state [ij]: the plaquette state is odd under 7r/2 rotations 
about the plaquette center. 

To measure the symmetry characteristic of the d-Mott 
state on a plaquette we propose an operator, Dp, which, 
when operating on a state, rotates a plaquette by 7r/2 
about its center or, alternatively, translates each par- 
ticle by one site along the four-site "ring". More con- 
cretely, for a particular basis state for demonstrative pur- 
poses, the single four-particle four-site plaquette (such 
as the lattice labeled in Fig. [Hb)) occupation could be 
|ni|n2fn3|n4-f ; ni|n2|n3|n4|) = |1100;0110) where sites 
1 and 2 are occupied with spin-up particles and 2 and 3 
are occupied with spin-down. The result of operating Dp 
on this state is Dp|1100;0110) = |0110;0011). The par- 
ticle at site 1 moves to 2, 2^3, 3^4, and 4^1 (the last 
through periodic boundary conditions on the plaquette). 
It is easy to see that the possible eigenvalues of Dp are 
±1, ±i since {Dp)"" ^ 1. 

For t' = Q and U/t ~ the ground state is six-fold 
degenerate and the operation oi Dp yields ±1 or ±i for 
each degenerate ground state. When t' = Q and U/t > 0, 
the degeneracy is broken and a gap A opens up. The 
operation of Dp onto the ground state now yields an 
eigenvalue equal to —1 indicating the existence of d- wave 
symmetry-the state has a non-zero value of the d-Mott 
operator {Dp). In Fig. [21 we plot the gap A with respect 
to U/t. The gap is zero for U = due to the degeneracy 
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FIG. 2: The gap A as a function of U/t for a single plaquette 
{t' — 0). The inset shows, approximately, the region in {U/t)- 
[t' /t) space where the d-Mott state is expected to exist. 



of the ground state while with finite U the degeneracy is 
lifted and the gap becomes nonzero. However, the gap 
also approaches zero in the limit U/t oo where there is 
no phase coherence between the atoms at different sites 
and the ground state is again degenerate. At a finite 
U/t « 6, the gap is maximum. 

We now consider what happens when t' ^ 0. An im- 
portant point to note is that Dp does not commute with 
the Hamiltonian for t' ^ 0, therefore eigenstates of H will 
not be eigenstates of Dp, altough H retains a global C4 
symmetry about the center of any one plaquette. How- 
ever, there are clearly defined regimes for non-zero values 
of t'. As mentioned, for t' — and U/t > 0, the ground 
state is gapped and it is an eigenstate of H and D p (with 
eigenvalue -1). For t' ^ A the ground state is approxi- 
mately a product state of single plaquette ground states 
and, hence, remains approximately a simultaneous eigen- 
state of H and Dp (by this we mean that the ground 
state will have the largest amplitude from basis states 
that are Dp eigenstates with eigenvalue —1). Thus, it is 
conceivable that the average value of Dp for any given 
plaquette will continue to be nonzero and close to — 1 
for a non-zero range of t' . On the other hand, when 
t' ^ t, the Hamiltonian is approximately the standard 
Hubbard model. Thus, the ground state is an essentially 
equal superposition oi Dp eigenstates, so one expects the 
value of (Dp), for any given plaquette, to be zero, since 
l + i-l) + i + i~i)=0. 

The question now becomes exactly what happens for 
U/t > and t' ^ A when the particles are first able 
to hop to the neighboring plaquettes. To answer this 
question we exactly diagonalize H for eight particles on 
an eight-site square lattice ladder with periodic boundary 
conditions in the "y" -direction depicted in Fig.[TJb) and 
vary both U/t and t' /t from zero to t' /t = 1 and beyond. 
(Note that a fully periodic system-torus-requires 16 sites 
with a quite large Hilbert space-~ 100 million states- 
which is not necessary for establishing the qualitative and 
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FIG. 3: (Color online) Ground state expectation value of the 
d-Mott operator, {Dp), as a function of t' /U (left) and t'/A 
(right) for the eight-site square lattice ladder system for the 
checkerboard Hubbard model at half filling for various values 
of U/t. 



semi-quantitative conclusions presented here.) The d- 
Mott character is examined by calculating the average 
{Dp). In Fig. [3] we plot {Dp) as a function of t' /U (left 
panel) and i'/A (right panel) for various values of U/t 
(weakly to strongly interacting). From the left panel we 
see that the behavior is largely universal when plotted 
versus t' /U and the value drops well below —1 for t' « 
t. However, the physics is better elucidated in the right 
panel. Here it is clearly seen that once t' becomes equal 
to, and exceeds, the single plaquette gap A the d-Mott 
order is quickly lost. Thus, the optimal parameter regime 
for observing the d-Mott state is at U/t w 6, where the 
gap A reaches the maximum and the d-Mott state, with a 
non-zero expectation value of the local plaquette operator 
Dp, is most stable against the inter-plaquette coupling 
t' . See the inset of Fig. [5] for a plot showing the region 
in {U /t)-{t' /t) space where the d-Mott state is expected. 

We now turn to how, in a cold atomic gas, the d-Mott 
state could be observed via time-of-flight measurements. 
First, we suggest a method to generate the d-Mott state 
in experiments. To begin, the potential depth of the 
short period optical lattice is ramped up to a large value 
Vs = 20Ej) to form a Mott state with one atom per lat- 
tice site [g], where Ep is the photon recoil energy. The 
potential depth of the long period optical lattice is then 
ramped up to V; = YIEp. In the Mott state, we have 
t — t' K, Q. To create the d-Mott state, the potential 
depth of the short period optical lattice is then ramped 
down to a small value (for instance, Vs = SEp) to en- 
hance the tunneling t inside the plaquettes. The poten- 
tial barrier between plaquettes is Vi + Vs — 20Ep, there- 
fore t' « 0. There are now four atoms in each plaquette, 
which corresponds to the half filled case discussed above. 
Non-zero t' can be obtained by lowering Vi and the inter- 
action strength U can be adjusted in experiments using 



Feshbach resonance [6|. With this method, various pa- 
rameter regimes can be reached with different U /t and 
t'/t. For the regime of U/t w 6, t'/A < 1 (i.e. t' /t < 0.3 
from Fig. [2]), a stable d-Mott state is obtained. We now 
address the question of the experimental observation of 
this state. 

Interestingly, the phase coherence in each plaquette in 
the d-Mott state can be experimentally identified in the 
time-of-flight (TOF) measurements. The d-Mott state 
in a single plaquette can be written as |x) = J2a fa |q^) 
where \a) is a basis state such that the action of Dp on 
a state rotates it to another basis state and the d-wave 
symmetry of the ground state wavefunction has been in- 
corporated into the eigen-coefficients /„. In the TOF 
experiment, we can measure the density of the spin up 
and spin down atoms separately. The TOF density dis- 
tribution for spin up atoms is 
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where Q (r) = mv/ht (m and t being the atomic mass 
and time of measurement, respectively), / is the index of 
plaquettes, 
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is the wavefunction of spin up atoms in plaquette I after 
TOF, b\ is the creation operator of the d-Mott state in 
the plaquette /, R/ is the position of the center of pla- 
quette /, j — 1,2,3,4 correspond to four lattice sites at 
each plaquette (see Fig. |4l^a)) with Si = {—a/2,— a/2), 
S2 = (a/2, -a/2), S3 = (a/2, a/2), S4 = (-a/2, a/2), a is 
the lattice spacing, fiaj is the spin state of atoms at site 
j and basis state a, and 5> (Q) is the Fourier transform of 
the Wannier function at each site. In the d-Mott state, 
we have (b\bii \ oc J//', yielding 
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In a single plaquette, we find that the d-Mott ground 
state can be written as 



\XP) - y(-i^p)*-M(A-7)|1100;0011) 

1=1 

+7IIIOO; 1100) - C{|1100; 1010) + |1100; 0101)}] 
2 

+ Y^{-Dpy-^[^{-\l{)l{);llOO) - |1010;0110) 



-|1010; 1001) -F|1010; 1010)} 
-2(A-7)|1010;0101)] 
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where A, 7 and ^ are parameters determined by {U, t, t') 
and satisfy normalization (16^^ -f 12 (A — 7) -I- 47^ — 1). 



FIG. 4: (Color online) (a) A single plaquette I. (b) Density 
distribution of d-Mott state in a time of flight experiment. 
The unit for Qx, and Qy is Tr/a. U/t = 5. 



When U oo, A \/3/6, ^ ^ 0, 7 ^ and the state 
becomes 

Ixp) ^(|1100; 0011)- 10110; 1001) + 10011; 1100) 



|1001;0110)) 
10101; 1010)) 



V3, 
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which has no double occupancy. Under the operation 
of Dp onto \xp) (either Eq. [5] or Eq. [6]), it is seen that 
an eigenvalue of — 1 is obtained. Using the wavefunction 
(Eq.O and Eq. ([3]), we find 



(n| (Q (r))) cx 4A^ |cos {QxO) — cos {Qya)\ 



(7) 



Note that the ground state of H (Eq. [T] with finite t') is 
approximately the d-Mott state when (Dp) « — 1 and, 
hence, it is approximately described by a product state 
of single plaquette ground states (Eq. [5]) with a TOP 
density distribution given by Eq. [71 

In Fig.Hfb), we plot the density distribution for the d- 
Mott state in TOF measurements with the image showing 
interference peaks at Qx + Qy = Nn/a, where N is an 
odd integer and Qx and Qy are multiples of ir/a. We 
want to emphasize that these peaks disappear in the TOF 
for a pure Mott state where one should observe a flat 
background, and this would be a way to contrast the 
present state from a regular Mott insulator. Further, the 
behavior of (Dp) as a function of t' (Fig. [3]) is exactly 
indicative of the "d-Mott"-ness of the ground state of 
H with a TOF signature shown in Fig.[4t^b). Therefore, 
monitoring the TOF interference peaks while varying U/t 
and t'/t would be an effective way to observe the d-Mott 
state. 

In Fig. [5l we plot the coefficient 4A^ with respect to 
the onsite energy U and find that the signal strength in- 
creases reaching 4/3 for U 00. However, the d-Mott 
state signal may not be observable for large U because the 
gap A is small at large U (see Fig. [5]), where a small inter- 
plaquette tunneling t' may mix the d-Mott state with 



FIG. 5: Signal intensity 4A'^ versus U /t. 



other states (see Fig. |31[b)), destroying the order and re- 
moving the interference signal. This is expected because 
for U ^ 00, the system is in a pure Mott state where no 
interference peaks are expected. Another possible road- 
block in the path of observing the interference pattern is 
that the temperature of the system must also be below 
A so that thermal fluctuations do not destroy the d-Mott 
order. But using U/t w5-6 and keeping fcsT ^ A should 
lead to the experimental observation of the d-Mott state. 

To conclude, we have described a scenario where a new 
state of matter, the d-Mott insulator, could be observed 
in a fermionic optical lattice via TOF measurements. 
Furthermore, we have provided semi-quantitative param- 
eter values where the state should most likely be seen 
{U/t « 6 and t' /t < 0.3). This could also serve as a rig- 
orous benchmark in further experimental studies toward 
the realization of the general two-dimensional square lat- 
tice Hubbard model in relation to its utility in the study 
of high-Tc cuprate superconductors. We believe that 
the direct TOF observation of the hitherto-unobserved 
d-Mott quantum insulator in a fermionic optical lattice 
would go a long way in establishing cold atomic gases as 
a useful optical lattice emulator of novel strongly corre- 
lated phases. 

This work is supported by ARO-DARPA. 
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